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We show that a heavy quark moving sufficiently fast through a quark-gluon plasma may lose energy 
by Cherenkov-radiating mesons. We demonstrate that this takes place in all strongly coupled, large- 
Nc. plasmas with a gravity dual. The energy loss is exactly calculable in these models despite being 
an C'(l/A'^c)- effect. We discuss phenomenological implications for heavy-ion collision experiments. 



1. Introduction. A remarkable conclusion from the 
Relativistic Heavy Ion Collider (RHIC) experiments [l[ is 
that the quark-gluon plasma does not behave as a weakly 
coupled gas of quarks and gluons, but rather as a strongly 
coupled fluid This makes the study of the plasma a 
challenging task. 

Experimentally, valuable information is obtained by 
analyzing the energy loss of energetic partons created 
in hard initial collisions. In order to use this information 
to learn about the plasma, a theoretical, quantitative un- 
derstanding of the different mechanisms of parton energy 
loss is needed. Several such mechanisms have been pre- 
viously studied, both in QCD itself [sl] and in the context 
of the gauge/gravity duality 4]. 

In this letter we will uncover a new mechanism whereby 
a sufficiently fast heavy quark traversing a strongly cou- 
pled plasma loses energy by Cherenkov-radiating in- 
medium mesons. We will first show that this takes place 
in all strongly coupled, large- A^,. theories with a gravity 
dual. Next we will calculate the energy loss in a simple 
example. Finally, we will discuss possible implications 
for heavy-ion collision experiments. 

2. Universality of the mechanism. This follows from 
two universal properties of the gauge/gravity duality (in 
the limit N^,X ^ oo): (i) the fact that the gauge the- 
ory deconfined phase is described by a black hole (BH) 
geometry and (ii) the fact that a finite number of 
quark flavours TV, is described by Ni D-brane probes [6| 
- see fig. [H In addition to the gauge theory directions, 
the gravity description always includes a radial direction 
which is dual to the gauge theory energy scale. The ra- 
dial position of the horizon is proportional to the plasma 
temperature T. The D-branes extend in the radial di- 
rection down to a minimum value proportional to the 
(constituent) quark mass M^. 

For sufficiently large M^/T, the D-branes sit outside 
the horizon [7H9|. In this phase, low-spin gauge theory 
mesons are described by small, normalizable fluctuations 
of scalar and vector flelds propagating on the branes, 
whose spectrum is discrete and gapped. In particular, 
this means that sufficiently heavy mesons survive decon- 
finement, in agreement with lattice and potential model 
predictions for real- world QCD [lo| . 
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FIG. 1: D-branes and open string in a BH geometry. 
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FIG. 2: DR for the transverse (black, continuous curve) and 
longitudinal (red, dashed curve) modes of a heavy vector me- 
son with Uiim = 0.35 in the D3/D7 system. The blue, contin- 
uous straight line corresponds %o lo = vq with < n < 1. 



Let w(q) be the in- medium dispersion relation (DR) 
for these mesons. As an illustrative example, the DR 
for vector mesons in the D3/D7 system is depicted in 
fig. m As q — > oo, the DR becomes linear: a;((7) ^ 
Wii^g, with Wii^ < 1. This subluminal limiting velocity, 
which is the same for all mesons, is easy to understand in 
the gravitational description ll| . Since highly energetic 
mesons are strongly attracted by the BH, their wave- 
function is very concentrated at the bottom of the branes. 
Consequently, their velocity is limited by the local speed 
of light at this point (see fig. [Ij. Because of the BH 
redshift, Wii„, is lower than the speed of light at infinity. In 
the gauge theory this translates into the statement that 
Wiin, is lower than the speed of light in the vacuum 12| . 

Consider now a heavy quark in the plasma. In the 
gravitational picture, this is described by a string that 
starts on the D-branes and falls through the horizon - 
see fig. [T] In order to model a highly energetic quark we 
consider a string whose endpoint moves with an arbitrary 
velocity v at an arbitrary radial position tq, where tq is 



2 



mesoi^^ 
quark _^^"quark 

FIG. 3: Effective quark-meson coupling. 



inversely proportional to the size of the gluon cloud that 
dresses the quark [14]. 

Two simple observations now lead to the effect that 
we are interested in. The first one is that the string end- 
point is charged under the scalar and vector fields on the 
branes. In the gauge theory, this corresponds to an effec- 
tive quark-meson coupling (see fig. 3) of order ^ l/\/Wl. 
The dynamics of the branes+string endpoint system is 
thus (a generalization of) that of classical electrodynam- 
ics in a medium in the presence of a fast-moving charge. 
The second observation is that the velocity of the quark 
may exceed the limiting velocity of the mesons, since the 
redshift at the position of the string endpoint is smaller 
than at the bottom of the branes. As in ordinary elec- 
trodynamics, if this happens then the string endpoint 
loses energy by Cherenkov-radiating into the fields on the 
brane. In the gauge theory, this translates into the quark 
losing energy by Cherenkov-radiating scalar and vector 
mesons. The rate of energy loss is set by the square of 
the coupling, and is therefore of order 
3. A quantitative example. In this section we will 
calculate the rate of energy loss in four-dimensional, 
SU{N,), Af = 4 super Yang-Mills (SYM) theory coupled 
to one quark flavour (but the result is valid for arbitrary 
Ni, see sec. 5). The dual description consists of a D7- 
brane probe in the supergravity background of black 
D3-branes. Following [121 we write the induced metric 
on the D7-brane worldvolume as ds^ = L^ds^{g) with 
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where = + , f = 1 - l/p^ / = 1 + l/p^ 
R = dR/dr, and — {t,x^} are the four gauge the- 
ory directions. R{r) describes the D7-brane embed- 
ding, with R{oo) — 2M^I \fxr . The dimensionless co- 
ordinates above are related to their dimensionful coun- 
terparts (denoted with tildes) through x^ = -kT i,^ , 
{r,i?, p} = {f,i?,/5}/7rL2r. 

The terms in the brane-t-string action relevant to our 
calculation are 
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where Fab ~ <9[a^6] and cr" — {a;^,r, fis}. The first term 
comes from expanding the Dirac-Born-Infcld part of the 
D7-brane action to quadratic order in the gauge field. 
The metric g that enters this term is that in eqn. (fTjl. 
which contains no factors of L; these have been absorbed 
in the definition of e in eqn. ([3]). The Wess-Zumino part 



of the D7-brane action will not contribute to our calcula- 
tion. The second term in ([2]) is the minimal coupling 
of the endpoint of an open string, whose worldline is 
parametrised by cr°(T), to the gauge field on the branes. 
We have omitted a similar coupling to the scalar fields, 
which will be considered in [15|. The coupling constant 
in (HI) is 
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(3) 



where Tdt = l/5s(27r)^^f is the D7-brane tension. As 
expected, e is of order l/\/N^, which justifies our neglect 
of terms of order higher than quadratic in the action. 

The second term in Q may be written as 
—e J d^aAaJ"". For simplicity, we will assume that the 
quark moves with constant velocity along a straight line 
at constant radial and angular positions, so we write 

r = 5(3) _ ^(3) _ X 

In reality, tq and v will of course decrease with time be- 
cause of the BH gravitational pull and the energy loss. 
However, for simplicity we will concentrate on the initial 
part of the trajectory (which is long provided the initial 
quark energy is larg e) for which ro and v are approxi- 
mately constant |l6j - see fig. [T] 

The rate of quark energy loss is given by minus the 
work per unit time done by the gauge field: 
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(5) 

Since real- world QCD has no internal S^, we focus on 
modes with no angular momentum on the S^. These take 
the form (a;'' , r) , (x" , r) , 0. We set = 

by a gauge choice. Further, we work with the Fourier- 
space components A^(aj,(j',r) and choose q = ((7,0,0), 
V = v{cos6,v sm9,0). After integrating over the S^, the 
relevant Fourier-space components of the current are 

= 27r6{uj-qvcos9)6{r^ro) x (1, cos6', sin6', 0) . (6) 

With this choice the only transverse mode of the gauge 
field excited by the source is A — A2. The equation of 
motion for this mode is 
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where J = 3^ , e ~ e/Vl-^ and Vl-^ ~ 2tt'^ is the volume of 
a unit S^. We solve ^ by expanding A as 



-4(w, g, r) = ^ X(w, q) £,n{q, r) 



(8) 



in terms of a basis of normalizable eigen-functions 
{£,n{q,r)} in the radial direction. These are solutions 
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of eqn. ([7]) with J = Q with g-dependent eigen-values 
Lj — ijjn{q), and satisfy the orthonormahty relations 
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Inserting the expansion ([8]) in ([7]), and using the eigen- 
state equation and the orthonormality relations, we find 



(10) 



where 



Jniuj,q) = J drj{uj,q)^niq,r) 

= 2TrS {uj — qv cos9) vsinO ^n{q,ro) ■ (11) 

Through the expansion (|5]) we have 'Kaluza-Klein' re- 
duced the five-dimensional gauge field to a discrete, in- 
finite tower of independent four-dimensional gauge fields 
{An{u}, q)}. Each of these fields is characterized by a q- 
dependent radial 'wave- function' r), as well as by a 
DR Lu = iOn{q), and couples to the quark with an effective 
strength erff(g,ro) = e£,n{q,ro). 

With retarded boundary conditions, as appropriate for 
the reaction to the quark's passage, eqn. ([TU]) yields 
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We now evaluate ([5]) to obtain the energy deposited 
on the n-th transverse mode. We first express 
^02(^1 wt, ro, r^o) as an integral over its Fourier compo- 
nents. We then integrate over frequencies trivially be- 
cause of the delta- function in pT|) . Finally, we set 
d^q — 2iTq^dqds, where s = cos 9, to arrive at 
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where s„ = Vn/v and v^^q) = uJn{l)/l is the phase veloc- 
ity of the mode. The Heaviside function confirms the ex- 
pected result: The quark only radiates into modes with 
phase velocity lower than v - those to the right of the 
dashed, vertical line in fig. [21 The numerical result for 
n = is plotted in fig. HI For fixed rg, the energy loss 
increases monotonically with v up to the maximum al- 
lowed value of w - the local speed of light at tq. As 
ro decreases, the characteristic momentum g<,i,ar of the 
modes contributing to the integral increases. As rg — >■ 
these modes become increasingly peaked at small r, and 
eeff(q'char, '"0) snd the energy loss diverge [11]. However, 
this mathematical divergence is removed by physical ef- 
fects we have not taken into account. For example, for 




FIG. 4: Energy loss into A(n=o) for an embedding with 
R{oo) = 1.32. The continuous curves correspond to ro = 
0.86, 0.97, 1.10, 1.25, 1.45. The dotted curve is defined by the 
endpoints of the constant-ro curves. 



sufficiently large q the radial profile of the mesons be- 
comes of order the string length and stringy effects be- 
come important [l3| . Also, mesons acquire widths F oc 
at large q [20|] and can no longer be treated as well defined 
quasiparticles. Finally, the approximation of a constant- 
V, constant-ro trajectory ceases to be valid whenever the 
energy loss rate becomes large. 

4. Phenomenology. The Cherenkov radiation of 
mesons by quarks depends only on the qualitative fea- 
tures of the DR of fig. m which are universal for all 
gauge theory plasmas with a dual gravity description [l3| . 
Moreover, it is conceivable that they may also hold for 
QCD mesons such as the J/i/' or the T (see e.g. the dis- 
cussion in 18]). Here we will examine some qualitative 
consequences of this assumption for HIC experiments. 

The energy lost into mesons would be reflected in a 
reduction of the heavy quark nuclear modiflcation factor 
21|. This would only occur for high enough quark 
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velocities, thus yielding a very particular behaviour of 
Raa- Note that the minimum quark velocity at which 
the reduction starts to occur may actually be higher than 
tiii„, since the quark energy must be larger than the mass 
of the radiated meson. For example, for a charm quark 
to radiate a J/i(j meson this condition yields v > 0.87. 
In fact, our calculation applies strictly only in the limit 
of infinite quark energy, which suggests that it should be 
more relevant to HIC experiments at the Large Hadron 
Collider than at RHIC. 

The radiated mesons would be preferentially emitted 
at a characteristic Cherenkov angle cos 0c = Vn^/v. Tak- 
ing the gravity result as guidance, v^m could be as low as 
~ 0.35 at the meson dissociation temperature [l^ . 
corresponding to an angle as large as 9^ « 1.21 rad. This 
emission pattern is similar to the emission of sound waves 
by an energetic parton 23[ in that both effects lead to a 
non-trivial angular structure. One important difference, 
however, is that the radiated heavy mesons would not 
thermalize and hence would not be part of a hydrody- 
namic shock wave. As in the Mach cone case, the meson 
emission pattern could be refiected in azimuthal dihadron 
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correlations triggered by a high-p^ hadron. Due to sur- 
face bias, the energetic parton in the triggered direction 
is hardly modified, while the one propagating in the op- 
posite direction moves through a significant amount of 
medium, emitting heavy mesons. Thus, under the above 
assumptions, the dihadron distribution with an associ- 
ated J/ip would have a ring- like structure peaked at an 
angle 9 k, n — 0^. 

5. Discussion. Cherenkov emission of gluons in the 
context of heavy ion collisions has been considered in 
[i^ . where the in-medium gluons are assumed to have 
space-like DR. Although some of the underlying physics 
is similar, the mechanism we have discussed is different in 
two respects. First, the radiated particles are colourless 
mesons, not gluons. Second, the gauge/gravity duality 
provides a large class of completely explicit and calcula- 
ble examples in which this mechanism is realized. 

We have focused on the transverse modes of the gauge 
field. Since the vector mesons are massive, there is a 
similar energy loss into the longitudinal modes 



paran, B. Fiol, A. Guijosa, C. Manuel, A. Paredes, 
L. Patifio and P. Townsend for discussions. We are sup- 
ported by a Marie Curie PIEF-GA-2008-220207 (JCS) 
and 2009SGR168, MEG FPA 2007-66665-G02 and GPAN 
GSD2007-00042 Gonsohder-Ingenio 2010 (DF, DM). 



151 



We calculated the energy deposited on the branes by 
the string endpoint. Since the branes sit outside the 
BH, this energy must stay on the branes (in the limit 
A'^,A oo). Because total conserved charges must 
agree, this energy is the same as the energy lost by the 
quark in the gauge theory. The Gherenkov angle 0^ is also 
the same, since it is determined by kinematics alone. In 
contrast, extracting unintegrated gauge theory quanti- 
ties {e.g. the differential power spectrum) would require 
computing the boundary stress-energy tensor as in [24 1. 

The validity of our results requires not just large 7V^, 
but also strong coupling. In this regime the holographic 
mesons behave as elementary excitations, as opposed to 
composite bound states, up to energies of order \/\T ^ 
T j7l Il2|. Despite the strong coupling requirement, our 
results might apply to asymptotically free theories as long 
as: (i) they are sufficiently strongly coupled at the scale 
T, so that the in-medium meson DR shares the qualita- 
tive features of that in fig. [21 and (ii) there is a non-zero 
quark-meson coupling in the medium. 

The energy loss is A^t-independent because (at leading 
order) the string endpoint couples directly to the gauge 
field on only one of the A'f D-branes. The \/N^ scal- 
ing of the energy loss does not necessarily imply that 
the analogous effect (if present) is small in A'^ = 3 real- 
world QGD, in particular at high quark velocities. Fur- 
thermore, its characteristic geometry and velocity depen- 
dence may make it easily identifiable. 

We close with a comment on the energy loss of heavy 
mesons. At A'^ — > oo, these mesons experience no drag 
12l [25[ . At finite A',. , pointlike heavy mesons experience a 
drag of order 0(1 /A^_^^) 26]. Gherenkov radiation implies 
an 0(1/A"„) drag for fast excited mesons describable as 
a long string with both endpoints on the D-branes p7t . 
since each endpoint may radiate as an individual quark. 
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